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Abstract 

Objective: Resting-state functional MRI (rsfMRI) has been increasingly used for understanding brain functional architecture. 
To date, nnost rsflVlRI studies have exploited blood oxygenation level-dependent (BOLD) contrast using gradient-echo (GE) 
echo planar imaging (EPI), which can suffer from image distortion and signal dropout due to magnetic susceptibility and 
inherent long echo time. In this study, the feasibility of passband balanced steady-state free precession (bSSFP) imaging for 
distortion-free and high-resolution rsfMRI was investigated. 

Methods: rsfMRI was performed in humans at 3 T and in rats at 7 T using bSSFP with short repetition time {TR = 4/2.5 ms 
respectively) in comparison with conventional GE-EPI. Resting-state networks (RSNs) were detected using independent 
component analysis. 

Results and Significance: RSNs derived from bSSFP images were shown to be spatially and spectrally comparable to those 
derived from GE-EPI images with considerable intra- and inter-subject reproducibility. High-resolution bSSFP images 
corresponded well to the anatomical images, with RSNs exquisitely co-localized to the gray matter. Furthermore, RSNs at 
areas of severe susceptibility such as human anterior prefrontal cortex and rat piriform cortex were proved accessible. These 
findings demonstrated for the first time that passband bSSFP approach can be a promising alternative to GE-EPI for rsfMRI. 
It offers distortion-free and high-resolution RSNs and is potentially suited for high field studies. 

Citation: Cheng JS, Gao PP, Zhou lY, Chan RW, Chan Q, et al. (2014) Resting-State fMRI Using Passband Balanced Steady-State Free Precession. PLoS ONE 9(3): 
e91075. doi:10.1371/journal.pone.0091075 

Editor: Yu-Feng Zang, Hangzhou Normal University, China 

Received August 22, 2013; Accepted February 9, 2014; Publislied March 12, 2014 

Copyriglit: © 2014 Cheng et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported in part by Hong Kong Research Grants Council (HKU7837/11M). The authors have declared that no additional external 
funding was received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: Queenie Chan has declared affiliation to a commercial company (Philips Healthcare). This does not alter the authors' adherence to all the 
PLOS ONE policies on sharing data and materials. 

* E-mail: *ewu@eee.hku.hk 



Introduction 

Resting-state functional MRI (rsfMRI) examines the temporal 
correlations in blood oxygenation level-dependent (BOLD) signal 
in the absence of stimulus or task. The two most common rsfMRI 
analysis methods are seed-based analysis and independent 
component analysis (ICA). While the seed-based analysis tests a 
specific hypothesis by cross-correlating voxel-wise rsfMRI data 
with the time course within a priori defined region of interest [1], 
ICA provides a hypothesis-free decomposition of the data into 
spatially maximally independent components [2]. Regardless of 
the technique, a consistent observation is that the coherent 
patterns of low frequency fluctuations (<0. 1 Hz) in BOLD signal 
reflect a series of functional networks, the so-caUed resting-state 
networks (RSNs), in both humans [1,3,4] and animals [5,6,7]. 
Extensive efforts in this field provide new insights into both brain 
functional organization [8,9] and functional reorganization during 
disease, aging and learning [10,11,12]. 

To date, most rsfMRI studies have been conducted with 
gradient-echo (GE) echo planar imaging (EPI), which provides a 



good compromise between spatial and temporal resolution. 
Although highly successful in most brain regions, this technique 
has limitations that reduce image quality in areas with high 
magnetic susceptibility. The requirement for long echo time (TE) 
to achieve BOLD sensitivity makes the technique susceptible to 
signal dropout [13], while the use of EPI for high speed can 
introduce severe image distortion [14]. Therefore, rsfMRI in brain 
regions near susceptibility boundaries, such as human anterior 
prefrontal cortex and rat piriform cortex, can be difficult. 
Additionally, high field improves the signal to noise ratio (SNR) 
and sensitivity of rsfMRI [15], but it also increases the field 
inhomogeneites and susceptibility effects that cause image 
distortion and signal dropout. To overcome these limitations, a 
variety of methods have been proposed, including parallel imaging 
[16], spiral-in/out acquisition [17], 3D z-shimming [18] and so 
on. However, each method has its own penalties such as decreased 
SNR, blurring with spiral readout or reduced temporal resolution. 
The other limiting factor in EPI approach is the significant T2*- 
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weighting during the long readout that can cause warping and loss 
of spatial resolution [19,20]. 

Balanced steady-state free precession (bSSFP) sequence is an 
imaging technique that uses rapid radiofrequency (RF) excitation 
pulses combined with fully balanced gradient pulses during each 
excitation repetition time (TR). Its short TR and readout duration, 
together with high SNR, have the potential for fast, distortion-free 
and high-resolution functional imaging. For a comprehensive 
review, see Miller [21]. The original bSSFP task-evoked flVIRI 
studies used the steep magnitude /phase transition in the bSSFP 
olf-resonance profile [22,23,24]. Despite its functional potential, 
the use of transition band bSSFP was limited by the need for 
multi-frequency acquisitions to locate the narrow range transition 
band. Later, it was demonstrated that functional contrast could 
also be achieved by using the relatively large flat portion of the 
bSSFP profile [2.5,26,27,28,29]. The contrast mechanism of this 
'passband bSSFP' is largely similar to that of spin-echo (SE). 
Recent years have seen compelling applications of the passband 
bSSFP in distortion-free and high-resolution functional imaging, 
such as layer-specific fMRI to hypoxic challenges in the mouse 
retinal [30] and high-fidelity auditory tonotopy mapping [31]. A 
preliminary study also reported odor-evoked activity in the human 
olfactory bulb, a small structure in high-susceptibility region [32]. 

In this study, the feasibility of bSSFP rsfMRI was investigated in 
humans at 3 T and in rats at 7 T. RSNs derived from bSSFP 
images using ICA were compared with those derived from GE- 
EPI images. Intra- and inter-subject reproducibility of bSSFP 
rsfMRI was also examined. Moreover, RSNs of high-resolution 
bSSFP rsfMRI were evaluated with regard to their anatomical 
correspondence. Finally, distortion-free bSSFP rsfMRI was 
demonstrated in himian anterior prefrontal cortex and rat piriform 
cortex. 

Methods 

Ethics Statement 

The human study was approved by the institutional review 
board of the University of Hong Kong. The participants were 
informed of the aims of our study before MRI examinations. FuU 
written informed consent was obtained from each participant. The 
animal study was approved by the Committee on the Use of Live 
Animals in Teaching and Research of the University of Hong 
Kong. All animal experiments were carried out in strict 
accordance with the recommendations in the Guidelines for the 
use of Experimental Animals of the University of Hong Kong. 

Human Imaging 

Eleven healthy right-handed subjects (6 male and 5 female; age 
23.3±2.8) participated in this study. MRI was performed on a 
Philips 3T MRI Achieva scanner (Philips Healthcare, Best, The 
Netherlands) with an 8-channel SENSE head coil. During rsfMRI, 
subjects were instructed to rest with their eyes open. bSSFP scan 
was performed with alternating RF pulse (P = n), number of slices 
(Nslices)=l, slice thickness = 5 mm, field of view (FOV) = 
240x240 mm2, matrix size = 64x64 (reconstructed to 128x128 
by zero-fining), TR/TE = 4/2 ms, number of signal averaged 
(NSA) = 1 2 yielding a temporal resolution = 3 s and number of 
dynamic scans (Ndyn) = 100. An optimal flip angle (FA) of 30o was 
estimated with Ti= 1820 ms and T2 = 99 ms [33] by cos(- 
FAopJ = (Ti/T2-l)/(Ti/T2+l) to provide the maximum flat 
passband region in the bSSFP olf-resonance profile [34] . Volume 
shimming was performed prior to each bSSFP acquisition to avoid 
the narrow transition band and minimize banding artifacts. For 
comparison, GE-EPI scan was performed with TR/TE = 3000/ 



30 ms, FA = 90°, SENSE factor = 2, Nslices = 33, shce thickness/ 
gap = 5/0 mm, FOV = 240x240 mm^, matrix size = 64x64 (re- 
constructed to 128x128 by zero-filling) and Ndyn= 100. 2D Tj- 
weighted (TiW) fast spin-echo (FSE) images were acrjuired for 
anatomical referencing, with TR/TE = 450/10 ms, FOV = 
240x240 mm^, matrix size = 240x240 and slice thickness = 
5 mm. 

During each session, a whole-brain GE-EPI (Nslices = 33) scan 
was first performed. Based on the orthogonal scout images, the 
GE-EPI scan was positioned parallel to the anterior commissure - 
posterior commissure plane with the center slice c()\'ering the 
posterior cingulate cortex, and other slices covering intraparietal 
cortex, motor cortex, visual cortex and superior temporal cortex. 
Subsequent bSSFP (Nslices = 1) scans were performed with each 
scan matching its GE-EPI counterpart slice location. To evaluate 
the bSSFP rsfMRI approach, four experiments were designed as 
follows. First, to examine the feasibility of bSSFP rsfMRI, five 
bSSFP scans were acquired from one subject (Subject 1) within 
one session, with each scan covering the posterior cingulate cortex 
(PCC), intraparietal cortex, motor cortex, visual cortex and 
superior temporal cortex respectively, where RSNs were most 
commonly reported [35]. Second, to investigate intra-subject 
reproducibility, four more sessions were performed on different 
days on Subject 1 with a bSSFP scan covering the PCC in each 
session. To examine the inter-subject reproducibility, one session 
was performed on each often other subjects (Subjects 2-1 1) with a 
bSSFP scan covering the visual cortex in each session. Third, to 
examine its potential for high-resolution rsfMRI, bSSFP scans with 
matrix size = 128x128 covering the PCC were also acquired from 
Subjects 2-5 in the above-mentioned sessions. Temporal resolu- 
tion of 3 s was maintained by SENSE acceleration (SENSE 
factor = 2) togetiier witii TR/TE = 4.8/2.4 ms and NSA = 15. 
Finally, to demonstrate the feasibility of bSSFP for distortion-free 
RSN detection, one session was performed on Subject 1 with a 
bSSFP scan covering the anterior prefrontal cortex (aPFC) near 
the air- tissue interface. 

Animal Imaging 

Nine normal male Sprague-Dawley rats (200-250 g) were used 
in this study. Animals were initially anaesthetized with 3% 
isoflurane and then mechanicaUy ventilated with 1.2-1.5% 
isoflurane via oral intubation. To minimize motion, animals were 
secured in the prone position on a plastic holder with the head 
fixed with a tooth bar and plastic screws in the ear canals. Animals 
were kept warm with circulating water while respiration rate, heart 
rate, oxygenation saturation and rectal temperature were contin- 
uously monitored and vital signs were maintained within normal 
physiological ranges [36,37,38]. 

MRI was performed on a 7 T Bruker scanner with a maximum 
gradient of 360 mT/m (70/16 PharmaScan, Bruker Biospin 
GmbH, Germany), a 72 mm birdcage transmit-only RF coil and 
an actively decoupled receive-only quadrature surface coil. Scout 
T2-weighted (T2W) images were first acquired in three planes with 
a rapid acquisition with relaxation enhancement (RARE) pulse 
sequence. For rsfMRI, bSSFP scans were acquired with alternat- 
ing RF pulse (P = Jt), FA= 19° as optimized using Ti = 1754 ms 
and T2 = 45 ms [29], Nslices =1, slice thickness = 1 mm, TR/ 
TE = 2.506/1.253 ms, matrix size = 64x64 and NSA= 9 resulting 
in a temporal resolution = 1 .5 s. Multiple bSSFP scans were 
performed on three animals with imaging planes covering primary 
visual cortex (VC) (Bregma - 7.2 mm) witii FOV = 32x32 mm^ 
on five animals covering primary somatosensory cortex (SC) 
(Bregma - 3.2 mm) with FOV = 32 x32 mm^ and on one animal 
covering piriform cortex (PC) (Bregma 3.3) with 
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FOV = 26x26 mml Each resting-state scan had Ndyn of 300, 
lasting 7 minutes 30 seconds. For comparison, one GE-EPI scan 
with TR/TE= 1500/18 ms, FA = 30°, matrix size = 64x64, 
Nslices=10 and slice thickness/gap = 1/0 mm was performed 
on each animal with one slice covering VC or SC. No resting-.state 
GE-EPI scan was performed to cover PC due to the severe image 
distortion and signal dropout. Prior to rsfMRI scans local 
shimming was performed with a FieldMap based procedure [39] 
to avoid the narrow transition band and minimize banding 
artifacts in bSSFP. For anatomical referencing, 2D T2W RARE 
images were acquired at the same slice locations with TR/ 
TE = 4200/36 ms, matrix size = 256 x256, echo train length = 4 
andNSA = 2. 

Data Preprocessing 

The corresponding slices in a GE-EPI scan that matched any 
bSSFP imaging plane were extracted and treated as a single-shce 
scan prior to pre-processing and rsIMRI analysis. The first five 
volumes of each resting-state scan were discarded to avoid possible 
non-equilibrium effect in bSSFP and GE-EPI dynamic scans. For 
the human data, the following preprocessing steps were carried out 
using FSL v3.3 (FMRIB, Oxford, UK): motion correction [40], 
removal of nonbrain structures from GE-EPI using BET [41] and 
those from bSSFP using manually drawn masks, spatial smoothing 
using Gaussian kernel with full width half-maximum (FWHM) 
5 mm and temporal band-pass filtering with 0.001—0.1 Hz cutoff. 
For the rat data, similar preprocessing steps were performed: intra- 
animal image realignment by 2D rigid-body transformation using 
AIR v5.2.5 (Roger Woods, UCLA, USA), spatial smoothing using 
Gaussian kernel with FWHM 0.5 mm and temporal band-pass 
FIR filtering with 0.005-0.1 Hz cutoff using MATLAB (Math- 
Work, Natick, MA). 

Resting-state fIVlRI Analysis 

Resting-state fMRI analysis was carried out using independent 
component analysis (ICA) on single-slice bSSFP and GE-EPI scans 
individually. ICA was implemented in MELODIC v3.1 [42] with 
15 independent components (ICs) for the human data and in 
GIFT v2.0d (http:/ /www.nitrc.org/projects/gift/) with 10 ICs for 
the rat data respectively. A mixing matrix was estimated in the 
ICA containing spatially ICs and corresponding time courses. The 
spatial IC was then converted to z-scores, which reflected the 
degree to which the time course of a given voxel correlated with 
that of the IC [43] . The histograms of IC z-scores were plotted and 
examined. The z-maps were thresholded based on a histogram 
mixture modeling at a probability (P-value) threshold >0.5 by 
using an alternative hypothesis-testing approach, placing equal loss 
to both false negatives and false positives (false positive rate = 0.5) 
[42,44,45]. While the thresholds differed in z-scores, they 
corresponded to the same probability. Statistical comparisons 
were carried out to examine the sensitivity of bSSFP for rsfMRI, 
the intra-subject and inter-subject reproducibility. The number of 
pixels with z-scores above the probabiUty threshold and the mean 
z-scores of these pixels were statistically compared between the 
bSSFP and GE-EPI derived RSNs using two-tailed paired t-test 
with p<0.05 considered as significant. AH the results were 
presented as means ± standard deviations. 

To assess low frequency (<0.1 Hz) power distribution, each 
time course was converted to power spectrum using Fast Fourier 
Transform (FFT) after detrending. The power spectrum was 
normalized by dividing its values by its maximum. 

RSNs were identified by examining each ICs spatial location 
and power spectrum, which should be locahzed in anatomically 
meaningful structures and dominated by low frequency [46]. 



Interpretations of RSNs were achieved with reference to previous 
ICA studies [47,48] in conjunction with atlases for human brain 
(Harvard-Oxford cortical structural atias, http://www.cma.mgh. 
harvard.edu/) and rat brain [49]. 

Results 

Human Data 

Figure 1 illustrates the seven RSNs derived from bSSFP images 
and their GE-EPI counterparts acquired from the same subject at 
identical slice locations, with the former showing a close spatial 
and spectral resemblance to the latter, as well as the typical, well- 
known RSNs demonstrated by numerous rsfMRI studies (For an 
overview see: Smith et al., 2009). Fig. lA shows the intraparietal 
sulcus, known to be part of the dorsal attention network [50]. 
Fig. IB illustrates the precentral gyrus, or the primary motor 
cortex (MC). Fig. IC shows the posterior cingulate cortex (PCC) 
and the lateral parietal cortex, posterior part of a regular default 
mode network. Figs. ID-F represent the visual networks spKt into 
occipital, medial and lateral visual cortices (VCs). Fig. IG shows 
the superior temporal cortex from the auditory network. The 
power spectra of all of the above RSNs, derived from both bSSFP 
and GE-EPI, clearly exhibited low frequency dominance. 

The intra-subject reproducibility of bSSFP rsfMRI is demon- 
strated in Figure 2. Spatial patterns of one RSN of the same 
subject examined in five sessions on different days were found 
similar. To examine the inter-subject reproducibility, resting-state 
bSSFP and GE-EPI scans were collected from ten subjects. 
Figure 3 shows that similar RSNs for visual cortex were 
consistently identified across all the ten subjects. Furthermore, z- 
scores of all RSNs derived from bSSFP and GE-EPI images were 
plotted into histograms, showing clear Gaussian distributions. 
Based on a histogram Gaussian/ Gamma mixture modeling, each 
z-map was thresholded at the same probability threshold (P>0.5). 
No significant difrer(;iice was detected in the number of pix(;ls and 
mean z-scores above thresholds between the RSNs derived from 
the two imaging methods (Table 1). 

High-resolution bSSFP images were also acquired at in-plane 
resolution of 1 .88 X 1 .88 mm^ (matrix size = 1 28 x 1 28, reconstruct- 
ed to 240x240 by zero-fiUing). As shown in Figure 4, high- 
resolution bSSFP (middle column) images corresponded well to 
the TiW FSE images (left column). Furthermore, without any co- 
registration the RSNs derived (right column) exquisitely co- 
localized to the gray matter segmented from TjW images using 
FAST [51]. 

Figure 5 illustrates the distortion-free RSN mapping capability 
of bSSFP. Tiu; anterior prefrontal cortex (aPFC) was imaged by 
TiW FSE, bSSFP and GE-EPI sequences with identical geometry. 
AH images were overlaid with edges of the TiW reference (Fig. 5A). 
At aPFC, image distortion and signal dropout were obvious in the 
GE-EPI image even with SENSE factor = 2 (as pointed by arrows 
in Fig. 5A). The bSSFP image, however, showed an excellent 
overlap with the TiW reference. A RSN representing aPFC was 
detected from bSSFP (Fig. 5B, left). In contrast, of all those- derived 
from GE-EPI two independent components were found located at 
aPFC. They were considered, however, unreliable because they 
bordered on the regions suffered from severe image distortion and 
signal dropout (Fig. 5B, middle and right). 

Animal Data 

bSSFP approach for rsfMRI was performed and compared with 
GE-EPI in isoflurane-anesthetized rats. As shown in Figure 6, 
robust RSNs were derived from bSSFP and GE-EPI images 
covering SC from all the five animals and covering VC from all 
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Figure 1. Resting-state networks (RSNs) of human derived from bSSFP and GE-EPI images. RSNs Identified representing dorsal attention, 
primary motor cortex, posterior part of default mode network, occipital visual cortex (VC), medial VC, lateral VC and auditory cortex (A-G) using 
independent component analysis. bSSFP and GE-EPI images were acquired from the same subject with identical temporal resolution, spatial 
geometry and resting-state paradigm. Power spectra, converted from time courses of corresponding independent components using Fast Fourier 
Transform, clearly indicate low frequency dominance. Z-maps were thresholded using histogram mixture modeling at an alternative hypothesis 
threshold P>0.5. 

doi:10.1371/journal.pone.0091075.g001 



the three animals. The spatial patterns of each RSN were 
consistent across animals and largely comparable between the two 
approaches. Moreover, the two spectra yielded by them were both 
low frequency dominated. High power in previously reported low 
frequency range (below 0.04 Hz) [52] was observed in most rats, 
with occasional appearance of peaks between 0.04 Hz and 0.1 Hz. 
It should also be noted that GE-EPI images were affected by 
susceptibility artifacts at lower brain areas close to the ear canal 
and airway (as pointed by arrows in Fig. 6), whereas bSSFP images 
were in good agreement with rat brain atlas [49] . 

Distortion-free RSN mapping capability of bSSFP was also 
investigated at the piriform cortex (PC) or the primary olfactory 
cortex, which is for one of the most important sensory functions in 
rats [53]. Figure 7A shows the typical T2W RARE, bSSFP and 
GE-EPI images at identical slice location of one animal, with PCs 
outlined in green based on the atlas. Image distortion and signal 
dropout at the inferior and lateral PCs were observed in GE-EPI 
image (as pointed by arrows in Fig. 7A, right). Since PCs are 
embedded with rich macrovasculature such as rostral rhinal veins 
[54], the image distortion and signal dropout were likely due to the 
susceptibility effects near large vessels in addition to brain tissue 
boundary. bSSFP images were, however, free of such susceptibility 
artifacts (Fig. 7A middle) and two unilateral RSNs for left and right 
PCs were identified (Fig. 7B). 

Discussion 

Feasibility of passband bSSFP rsfMRI 

As demonstrated in this study, passband bSSFP sequence is a 
feasible alternative to GE-EPI for rsfMRI. As shown in Figs. 1,2,3 
and 6 and Table 1, RSNs derived from bSSFP and GE-EPI 
images were essentially similar, both in close resemblance to the 
typical, well-known RSNs reported by numerous rsfMRI studies 



[47,55]. Considerable intra- and inter-subject reproducibility was 
also observed (Figs. 2, 3 and Table 1), indicating the robustness of 
the proposed bSSFP approach for rsfMRI. 

The resemblance in RSNs derived from bSSFP and GE-EPI 
images likely stems from the BOLD sensitivity common to both 
methods. Comparing to GE-EPI, the exact contrast mechanism of 
bSSFP remains elusive. It depends on the protocol parameters, 
particularly TR and FA [56,57]. By employing large FA, bSSFP is 
optimized for the passband; while by employing small FA, bSSFP 
is optimized for the transition band. At relatively large FA (30°/ 
19°) and short TR (4.0/2.5 ms) used in the current human and rat 
studies respectively, the off-resonance effects around microvascu- 
lature created by deoxyhemoglobin are in the same spatial scale of 
water diffusion distance and therefore the dynamic averaging 
effect cannot be refocused, leading to signal increase during local 
neuronal activation. This mechanism can be considered similar to 
that of TgW spin echo (SE). The main difference from SE is that 
bSSFP sequence allows for stimulated echo pathways that can give 
rise to increase in signal and BOLD contrast. bSSFP signal may 
also include an intravascular component of macrovasculature, 
mainly due to water molecules diffusion in and out of the red 
blood cells [25]. bSSFP with lengthened TR for high spatial 
resolution exhibits Tj* dephasing and the sequence behaves 
essentially like a poorly-spoiled GE [58]. Therefore, the T2*W also 
contributed to the BOLD sensitivity in the present study of high- 
resolution bSSFP rsfMRI. 

The above bSSFP contrast mechanism was established in the 
context of task-evoked responses. Although much is known about 
the neurophysiology underlying task-evoked BOLD responses 
[59], the neurophysiological basis of spontaneous modulation of 
BOLD signal remains to be debated. For example, nitric oxide 
synthase inhibition blocks task-evoked blood flow increases, but 
spontaneous fluctuations in blood flow can remain unaffected [60] . 
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Figure 2. Intra-subject reproducibility of bSSFP for rsfiVIRI. RSN representing default mode networl< was consistently identified from bSSFP 
and GE-EPI scans acquired from one subject at five sessions. Z-scores were plotted into histograms, showing clear Gaussian distributions. Z-maps 
were thresholded using histogram mixture modeling at an alternative hypothesis threshold P>0.5. 
doi:1 0.1 371 /journal.pone.0091 075.g002 



Nevertheles.s, it has been shown that both spontaneous and task- 
evoked GE-EPI BOLD responses depend on T2* BOLD contrast 
[61]. Hence, the spontaneous modulation of bSSFP signal most 
likely depends on T2 contrast, the same as that of bSSFP task- 
evoked responses. In addition, the power spectra of RSNs derived 
from both bSSFP and GE-EPI were shown to be dominated by 
low frequency, indicating that there exists a universal frequency 
signature of the spontaneous modulation of BOLD signal at 
resting-state. 



Sensitivity of passband bSSFP rsfMRI is expected to be 
dependent of acquisition parameters. Parameters such as TR, 
TE and FA can influence the contrast mechanism by which the 
sensitivity changes. Changes in partial volume effects at different 
spatial resolution are also expected to impact sensitivity. However, 
when selecting the acquisition parameters the significance of the 
above effects has to be weighed against other constraints such as 
the flatness of passband, the width of passband (typically 
0.75 xTR ' Hz) and RF specific absorption rate. In the current 
study with very short TR and optimized large FA, no significant 
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Figure 3. Inter-subject reproducibility of bSSFP for rsfMRI. RSN representing visual cortex was consistently identified from bSSFP and GE-EPI 
scans acquired from ten subjects. Z-scores were plotted into histograms, showing clear Gaussian distributions. Z-maps were thresholded using 
histogram mixture modeling at an alternative hypothesis threshold P>0.5. 
doi:1 0.1 371 /journal.pone.0091 075.g003 



difference was detected in the number of pixels and mean z-scores 
above thresholds between the RSNs derived from bSSFP and GE- 
EPI data. Moreover, at submilhmeter resolution most RSNs 
derived from bSSFP exhibited more focal clusters than those from 



Table 1. bSSFP and GE-EPI rsfMRI reproducibility. 







bSSFP 


GE-EPI 


P 


Intra-subject reproduc 


bility (1 subject with 5 


sessions) 




Mean z 


3.79±0.35 


3.85 ±0.08 


0.30 


Number of pixels 


530±n9 


560 ±82 


0.75 


Inter-subject reproduc 


bility (10 subjects with 


single session) 




Mean z 


3.73±0.59 


3.97 ±0.48 


0.28 


Number of pixels 


585±161 


475±111 


0.10 



Means and standard deviations of the number of pixels with z-scores above the 
probability threshold {P>0.5) and the mean z-scores of these pixels in RSNs in 
Figures 2 and 3. Statistical analysis was performed using two-tailed paired t-test. 
doi:1 0.1 371/journal.pone.0091 a75.ta01 



GE-EPI with rapidly varying z-scores among adjacent pixels 
(Fig. 6). Since bSSFP functional contrast is composed of not only 
the extravascular component sensitive to microvasculature but also 
the undesirable intravascular component sensitive to macrovascu- 
lature [57], the focal clusters in bSSFP derived RSNs maps might 
be due to partial volume effects from large veins. To minimize 
such undesirable effects, slight diffusion weighting or inversion 
recovery can be combined with bSSFP acquisition [30]. Never- 
theless, it remains imperative to systematically investigate the 
sensitivity of passband bSSFP for rsfMRI in the future studies. 

High-resolution and Distortion-free RSN mapping 

In the presence of signal dropout and image distortion, even 
with high encoding resolution, the effective spatial resolution can 
still be greatly reduced in GE-EPI with neighboring voxels 
inaccurately assigned in the image space. Furthermore, high 
acquisition resolution lengthens the readout and can further distort 
the image. While In passband bSSFP, the rapid excitation 
repetition intervals allow short segmented readouts, resulting in 
images with distortion levels similar to that of anatomical 
acquisitions (see human aPFC and rat PC images acquired with 
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Figure 4. High-resolution RSNs derived from bSSFP images. T,W FSE anatomical references, bSSFP images at high-resolution (in-plane 
resolution 1.88x1.88 mm^ matrix size = 128x128, reconstructed to 240x240 by zero-filling) and zoomed and segmented brain images. High- 
resolution bSSFP images well corresponded to their T, W references, with RSNs exquisitely co-localized to the gray matter without any co-registration. 
T,W images were segmented into white matter, gray matter and cerebrospinal fluid. 
doi:l 0.1 371/journal.pone.0091 075.g004 



GE-EPI and bSSFP in Figs. 5 and 7A). Tlierefore, RSNs derived 
from bSSFP can be interpreted for functional connectivity with 
high fidelity. As sliown in Fig. 4, RSNs derived from high- 
resolution bSSFP (1.88x1.88 mm^) exquisitely co-localized to the 
gray matter without any co-registration (Fig. 4). bSSFP also 
detected a RSN at aPFC (Fig. 5B), which has not been reported by 
GE-EPI at this slice location likely due to image distortion and 
signal dropout. Note that such a pattern has been detected by a 
recent SE-EPI study at 7 T [62]. Nevertheless, bSSFP approach 
offers advantages over SE-EPI approach for rsfMRI. It yields less 
image distortion with reduced readout, and it may also provide 
higher sensitivity, with higher signal and BOLD contrast arising 
from the stimulated echo pathways [27,29]. In our rat results at 7 
T, two unilateral components localized in PC were derived from 
bSSFP (Fig. 7). Bilateral PC RSNs have been reported in previous 
rat studies using GE-EPI, in which image distortion and signal 
dropout were alleviated by lower field strength at 4.7 T [5] or 
parallel imaging technique [6]. The absence of unilateral 
components in those studies may be caused by ICA underestima- 
tion and/ or a high degree of temporal correlation between intra- 
and interhemispheric network. Homologous regions have been 
known to operate both uni- and bilaterally [63]. Note that 
unilateral components representing other sensory networks have 



been reported in somatosensory and motor cortices in rats and 
mice [6,47]. 

While various methods have been proposed to reduce signal 
dropout and image distortion with conventional GE-EPI 
[16,17,18,64], most of such approaches can be adopted for bSSFP 
imaging. Moreover, the inherently distortion-free nature of bSSFP 
circumvents the additional processes in both acquisition and post- 
processing that are often employed to tackle the image distortion 
problem. 

Technical Considerations 

First, bSSFP requires good shimming to avoid banding artifacts 
particularly at high field. In this study, localized shimming was 
carefully performed within the single image slab prior to resting- 
state acquisition to minimize the banding artifacts [29]. Second, 
this study demonstrated the bSSFP rsfMRI with single slice 
acquisition. The differences in RSNs (medial VC for Subjects 2—9, 
lateral VC for Subject 10 and occipital VC for Subject 11 shown 
in Figure 3) could be caused by slightly different shce locations 
besides inter-subject variations in functional connectivity. Howev- 
er, fast 3D volumetric bSSFP fMRI has been demonstrated using 
various 3D k-space trajectories [22,65] and with SENSE 
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TjW bSSFP GE-EPI 




B Independent Components 




Figure 5. A distortion-free RSN derived from bSSFP images at 
human anterior prefrontal cortex (aPFC). A: T,W FSE, bSSFP and 
GE-EPI images at identical slice location were overlaid by edges of T,W 
images detected using FSL-slicer tool. Image distortion and signal 
dropout in GE-EPI (SENSE factor = 2) were pointed by an arrow head, 
while bSSFP image showed excellent overlap with TiW images. B: 
Independent components derived from bSSFP (left) and GE-EPI (middle 
and right) situated at aPFC. Note that the clusters within the latter two 
were adjacent to image distortion and signal dropout, affecting the 
precise interpretation of their neurophysiological significance. 
doi:10.1371/journal.pone.0091075.g005 

acceleration [66] . Such 3D acquisition can be readily adopted for 
rsfMRI. To cover the whole brain, a previous bSSFP fMRI study 
combined two acquisitions with different phase cycling angles to 
remove the banding artifacts [27]. Yet, recent fMRI studies have 
simulated [67] and implemented [68] a frequency-modulated 
bSSFP technique [69] enabling dynamic acquisition with different 
phase cycling angles without reestablishing the steady-state. This 
frequency-modulated technique could enable the bSSFP rsfMRI 
to measure temporal correlations within full brain. Combining 
such frequency-modulated technique, parallel imaging and/or 
compressed sensing and 3D acquisition, high spatial and temporal 
resolution bSSFP with full brain coverage can be achieved for 
rsfMRI. 



T,W bSSFP GE-EPI 




B Left Piriform Cortex Right Pirifonn Cortex 




Figure 7. A distortion-free RSN derived from bSSFP images at 
rat piriform cortex (PC). A: TjW RARE, bSSFP and GE-EPI images at 
identical slice location covering PC. PCs were overlaid by green 
contours based on atlas [49]. Note that GE-EPI suffered from image 
distortion and signal dropout at inferior and lateral PCs as pointed by 
arrow heads, while bSSFP was in good agreement with T2W image. B: 
Unilateral RSNs for left and right PCs derived from bSSFP. PC, Bregma 
3.3 mm. 

doi:10.1371/journal.pone.0091075.g007 

Conclusion 

This study shows that passband bSSFP is feasible for rsfMRI 
studies free of severe image distortion and signal dropout. RSNs 
derived from bSSFP images were shown spatially and spectrally 
comparable to those derived from GE-EPI images with consider- 
able intra- and inter-subject reproducibility in humans at 3 T and 
in rats at 7 T. High-resolution bSSFP corresponded well to the 
anatomical images, with RSNs exquisitely co-localized to gray 
matter. Furthermore, RSNs at areas of severe susceptibility were 
proved accessible including human anterior prefrontal cortex and 
rat piriform cortex. These findings demonstrated for the first time 
that passband bSSFP approach can be a promising alternative to 
GE-EPI for rsfMRI. It offers distortion-free and high-resolution 
RSNs and is potentially suited for high field studies. 



Primary Somantosensory Cortex Primary Visual Cortex 

Animal 1 Animal 2 Animal 3 Animal 4 Animal 5 Animal 6 Animal 7 Animal 8 




Frequency (Hz) 



Figure 6. RSNs of rats derived from bSSFP and GE-EPI images. Consistent RSNs including primary somatosensory cortex (SC) and primary 
visual cortex (VC) across animals were obtained from both bSSFP and GE-EPI images with individual power spectra dominated by very low frequency 
(<0.04 Hz). SC, Bregma - 3.2 mm; VC, Bregma - 7.2 mm. 
doi:1 0.1 371/journal.pone.0091 075.g006 
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